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We describe a direct method to reconstruct the transverse proton source formed in a rel- 
ativistic heavy ion collision, making use of experimentally measured proton and deuteron 
^ ', spectra and assuming that deuterons are formed via two-nucleon coalescence. We show 
O I that an ambiguity with respect to the source temperature still persists and we indicate a 

' possible solution to the problem. 
0^ ■ 

^ I It has been recently shown that a simple description of the proton phase space dis- 
tribution can provide a good qualitative understanding of deuteron spectra [0]. On the 
other hand, it appeared clear that a more solid method to extract the properties of the 
source had to be established. We have therefore reconstructed the phase space distribu- 
. tion of protons directly from the observed proton and deuteron spectra, exploiting the 
coalescence prescription together with the notion of collective flow. Here we will outline 
k> I the procedure while all details can be found in 0. We make use of a relativistic descrip- 
■ tion of collective flow based on the boost-invariant picture for the longitudinal expansion, 
together with a longitudinally-independent transverse velocity, and therefore write the 
proton phase space distribution as 

fpix,p) = i27rfexpi-pX{x)/To)Bpnp{r^), (1) 

where p^u^{x) = 7(r) {m± cosh(?/ — rj) — p± ■ v{f)) is the energy in the global frame and 
Bp is the normalisation coefficient of the Boltzmann distribution in the local frame. The 
local density np(r) is assumed to be independent of the longitudinal rapidity. 

The deuteron phase space distribution is calculated using the coalescence model. Its 
evaluation is simplified when considering large and hot systems, neglecting the smearing 
effect of the deuteron Wigner density in comparison to the characteristic scales of the 
system in position and momentum space. One therefore obtains the deuteron phase space 
distribution 

/rf(a;,p)^|i?„,[/,(x,p/2)]\ (2) 

The neutron to proton ratio in the source was taken to be Rnp = 1.2. The deuteron phase 
space distribution has the same structure as the proton one in eq. (|I]), now with np{r) 
replaced by nrf(r) = Xdn^ir), where = 3/8 Rnp{27rf B'^/ Bd- 
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One can now calculate proton and deuteron invariant momentum spectra using the 
Cooper-Frye formula, on an approximate freeze-out hypersurface of constant longitudinal 
proper time tq. The integrations over the space-time rapidity and the azimuthal angle can 
be easily performed, yielding an expression in terms of the two functions v{r) and n{r). 
The ambiguity in the description of the single particle spectrum is explicit, since the two 
functions cannot be mapped out uniquely from only one function as the transverse momen- 
tum spectrum. To partially remove the ambiguity, we first change the integration variable 
and introduce the auxiliary function h through the relation vdvh{v) = r dr Ton{r), ob- 
taining the new expression for the momentum spectrum 

S{p^) =AnBm^ £dvv K, (^) Jo (^) h{v) . (3) 

The one-to-one correspondence between n{v) and S{p±) is now evident. We then make 
use of the coalescence model and the definition of n, both for protons and for deuterons, 
to obtain a first order differential equation that can be directly integrated and gives the 
closed solution 

.^ = 2^/*„.3M. (4) 

To nd[u) 

Therefore, by independently extracting the functions fip and from the observed mo- 
mentum spectra Sp{p±) and Sd{p±), we can find the function r{v) by a simple numerical 
integration. Inverting the obtained function as r{v) v{r), we obtain the collective 
velocity profile. We also obtain the local proton density as 

n {r) = ^ ^^^^^ ^P^^^^)) (5) 
^ r dr To 

The described procedure was applied to the transverse momentum spectra, resulting 
from Pb+Pb collisions and measured by the NA44 collaboration at the CERN-SPS [^]. 
We fitted hp and to these data, using eq. @ for protons and deuterons and assuming 
different values for Tq, from 50 MeV to 150 MeV. We chose a form of the profile functions 
characterized by three parameters, and we extracted their value from the experimental 
data with a Monte Carlo search minimising x^- The fitted spectra are shown, for the 
extreme values of temperature considered, in the top part of Fig.[l|. Although evaluated 
at different temperatures, they are indistinguishable from one another. On the other 
hand, the profiles h result to be very different for different temperatures. 

In all calculations the freeze-out time was fixed at tq = 10 fm/c After numerical 
integration of eq. (^, we obtained the function v{r) shown in Fig.|^ for the two extreme 
temperatures. It shows a linear rise at small r and saturates for large r. The velocity 
profiles clearly depend on the temperature chosen. The local proton density is also plotted 
in Fig.|^. At high temperature the density shows a shell-like structure which disappears 
as the temperature is lower. Similar shell-like structures have been recently found in an- 
alytic solutions of non-relativistic hydrodynamics . From the plots one can observe the 
different transverse sizes corresponding to different temperatures. It is therefore necessary 
to know To precisely in order to determine the system size. This information cannot be 
extracted solely from proton and deuteron spectra. Contrary to what is commonly done, 
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Figure 1. Comparison of the transverse spec- 
tra of p, d and t, calculated within the im- 
proved coalescence model, with the measure- 
ments obtained by the NA44 experiment [Q]. 



Figure 2. Reconstructed local density and 
flow of protons as function of the trans- 
verse radius, for the extreme values of Tq. 
Up is rescaled to be compared with v. 



source radii cannot be extracted from the d/p^ ratio {B2 ~ emitting volume) unless Tq is 
known. 

It is clear that to resolve the remaining ambiguity in the source temperature one needs 
some additional experimental information. Recently, tttt correlation data were used as 
a constraint Since pions may freeze-out in a different way than protons, it would 
be even better to consider pp correlations, although they are more sensitive to final 
state interactions than pions. On the other hand, heavier clusters can provide additional 
constraints. We now address this issue following and we describe the fusion process of 
A-nucleons into a bound state within the density matrix formalism. Making use of the 
same approximation leading to eq. (0), one obtains cluster phase space distribution 



\N 



fc{x,p)^9A'-^^[U{x,p/A)]\ (6) 

For N = Z = 1, A = 2 and (72 = 3/8 it reduces to eq. (|). The statistical prefactor is 
of crucial importance. In the conventional approach we have qa = (25*^ + 1)2'! A^! / 2^A\, 
which results in ga = 3/8, gt = gme = 1/12, and giue = 1/96. Although straightforward, 
this approach has been successful only for the description of deuteron production. In fact, 
it significantly underestimates the yields of heavier clusters. A possible improvement can 
be achieved by allowing for additional formation processes. Besides the already considered 
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direct process, it is possible that deuteron-like correlations contribute to cluster formation. 
For the triton, as an example, we should therefore account for the possibility that a proton 
and a neutron are already in a bound state with deuteron quantum numbers and coalesce 
with another neutron, with a different statistical prefactor. More precisely we can write 
Qt = gpnn^t + '^gpn^dddn^t, wherc the factor 2 counts the different ways to associate 
the proton with the two neutrons in forming a deuteron. The spin-isospin counting is 
straightforward and gives the modified statistical prefactor gt = 1/3, therefore enhancing 
the triton yield by a factor 4. The same arguments apply to ^He, so that ga^^ = 1/3. The 
case of^He is more involved. Counting all the possible processes we obtained g4He = 13/48, 
so that the ^He yield is increased by a factor 26. For more details see pf. 

We now use the coalescence model in this improved version to calculate the trans- 
verse mass spectrum of tritons. Using the flow and density profiles extracted from the 
analysis of p and d spectra, we examine to what extent heavier clusters can constrain 
the ambiguity of the temperature. Using eqs. (|l]) and (H), we obtain nt{r) = XtUp^r), 
with At = gt^RnpY /2 (271)^ Bp/ Bt, while the collective velocity is the same for all clus- 
ters. The triton spectrum is plotted in Fig. |l|, together with the previously fitted p and 
d spectra. The absolute values and the shape compare very well with the experimen- 
tal data. This confirms, that the improved statistical approach is consistent with the 
measured spectra. Furthermore, the results obtained with the two extreme temperatures 
show a different curvature. The high temperature case presents a clear bending over, 
absent for low temperature. We argue that this specific difference might narrow down 
the allowed temperatures and therefore provide additional constraints. Unfortunately, 
the triton spectrum was measured only in a limited range in transverse momentum, and 
therefore a quantitative fit is not useful. 

We also suggest that a further possibility to constrain the temperature lies in the com- 
bined study of single and composite spectra, together with pp correlations. The extracted 
profiles shown in Fig. |^ could in fact be used to evaluate the correlation function, which 
is sensitive to temperature and flow in a different way respect to inverse slopes of spec- 
tra. These remarks are important for the interpretation of the future experiments at the 
Relativistic Heavy Ion Collider. 
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